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i465/rac^-  Scintillations  in  the  intensity  of  an  acoustic  signal 
are  a  common  feature  of  propagation  of  sound  in  the  sea, 
manifesting  temporal  variability  in  the  index  of  refraction 
(sound  velocity)  of  the  medium.  In  this  paper,  a  recently 
developed  high-order  perturbation  method  is  described  and 
applied  to  the  problem  of  sound  propagation  in  the  sea.  The 
method  uses  a  “canonical”  solution  (sound  speed  profile)  to 
form  a  set  of  basis  functions  that  span  the  solution  space  and 
adequately  represent  the  exact  eigenvalue  problem.  The  basis 
functions  used  in  the  calculations  are  derived  from  sound 
speed  profiles  obtained  in  an  acoustic  propagation  experiment 
conducted  in  a  shallow-water  region  of  the  Mediterranean.  At 
particular  source  frequencies,  calculations  of  modal  functions 
and  acoustic  transmission  loss  were  compared  for  the  mean 
and  several  perturbed  profiles.  The  results  confirm  the  signifi¬ 
cant  effects  on  acoustic  transmission  of  seemingly  minor 
variations  in  sound  speed  and,  moreover,  demonstrate  the 
efficacy  of  the  new  perturbation  method  in  handling  such 
problems. 

INTRODUCTION 

The  propagation  of  sound  in  the  ocean  is  inevitably 
accompanied  by  fluctuations  in  the  amplitude  and  phase  of  an 
acoustic  signal  received  at  large  distances  from  the  source. 
The  fluctuations,  or  scintillations,  are  manifestations  not  only 
of  changing  patterns  of  interaction  with  the  bottom  and  sur¬ 
face,  particularly  important  in  shallow- water  propagation,  but 
also  passages  of  the  wave  through  time- vary  ing  inhomogene¬ 
ities  in  the  ocean  medium.  The  variability  in  acoustic 
propagation  can  be  considered  to  arise  from  variations  in  the 
index  of  refraction,  or  sound  velocity,  of  the  medium,  which, 
in  turn,  are  induced  by  a  variety  of  ocean  processes  covering 
a  wide  range  of  temporal  and  spatial  scales.  The  fluctuations 
in  the  acoustic  field  arc  more  than  “merely”  of  academic 
interest,  since  they  result  in  significant  degradations  in  the 
performance  of  underwater  acoustic  systems. 

Theoretical  approaches  to  the  propagation  of  sound  in  a 
realistic  ocean  environment  are  still  evolving.  Well-known 
deterministic  numerical  models  include  the  parabolic  equa¬ 
tion  (PE)  method,  fast  field  programs  (FFP),  the  normal  mode 
method,  finite  element  methods,  finite  difference  approaches, 
etc.  Attempts  to  understand  the  non-deterministic  component 
of  acoustic  intensity  fluctuations  have  been  centered  around 
the  theory  of  propagation  through  a  random  medium.  Because 
the  pertinent  stochastic  wave  equation  cannot  be  solved 


deterministically ,  solutions  are  obtained  in  terms  of  statistical 
moments  of  the  acoustic  pressure,  such  as  the  expected  value 
(ensemble  average),  various  correlation  functions,  and  higher 
order  moments.  The  fourth  moment  is  of  particular  interest 
since  this  gives  the  intensity  fluctuations.  Much  of  the  new 
work  in  extending  the  range  of  validity  to  large  intensity 
fluctuations  arising  from  multiple  scattering  is  based  on  the 
non-deterministic  version  of  the  PE  method.  The  PE  method 
is  a  simplification  of  the  wave  equation  for  the  case  of 
scattering  that  is  mostly  in  the  forward  direction.  The  two  main 
methods  of  solving  the  non-deterministic  PE  are  the  Feynman 
path  integral  (used  by  Flatte,  et  al  [1])  and  the  moment- 
equation  (used  by  Uscinski,  eial  [2]).  Both  methods  have  had 
some  recent  success  in  attributing  measured  intensity  fluctua¬ 
tions  to  the  action  of  internal  waves  in  deep  water.  The 
preceding  success  has  been  possible  because  of  the  availability 
of  the  required  statistical  description  of  the  stochastic  me¬ 
dium,  viz,  the  Garrett-Munk  (GM)  empirical  model  [3,4].  In 
the  upper  ocean  and  in  shallow  water  the  GM  model  turns  out 
to  be  inappropriate;  instead,  a  quasi-deterministic  explanation 
involving  soliton  propagation  appears  to  be  more  relevant  [5]. 

In  this  paper  a  recently  developed  high  order  perturbation 
technique  is  applied  to  the  problem  of  sound  propagation  in  a 
varying  underwater  acoustics  environment.  The  method  uses 
a  “canonical”  solution  (sound  speed  profile)  to  form  a  set  of 
basis  functions  that  span  the  solution  space  and  adequately 
represent  the  exact  eigen  value  problem.  This  new  “perturbation 
method”  is  far  less  limiting  than  conventional  pertur¬ 
bation  approaches  since  higher  order  terms  are  not  neglected. 
The  basis  functions  used  in  the  calculations  are  derived  from 
sound  speed  profiles  obtained  in  an  acoustic  propagation 
experiment  conducted  in  a  shallow-water  region  of  the 
Mediterranean.  The  question  as  to  whether  the  varying  sound 
speed  profiles  are  manifestations  of  stochastic  phenomena  or 
quasi-deterministic  events  in  the  ocean  environment  is 
secondary  to  the  purpose  at  hand.  Simply  stated,  our  primary 
objective  is  to  demonstrate  that  the  perturbation  technique 
developed  is  eminently  suitable  for  handling  problems  such  as 
those  treated  here  and,  in  fact,  offers  considerable  advantages 
over  conventional  perturbation  techniques. 

The  paper  is  organized  as  follows.  In  the  following  section 
the  newly  developed  perturbation  theory  is  summarized.  This 
is  followed  by  a  brief  description  of  the  measured  sound  speed 
profiles  and  the  environment  in  which  they  were  obtained. 
Next,  examples  of  the  computed  acoustic  transmission  loss  for 
several  mean  and  “perturbed  profiles”  arc  discussed.  Finally, 
conclusions  are  provided. 
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2.  THEORY 

2.1  A  New  Full-Perturbation  Method  for  Normal  Modes 

In  this  section  we  outline  a  new  full  perturbation  approach 
to  obtain  solutions  to  the  normal  mode  problem.  This  method 
proves  to  be  very  fast  and  particularly  suitable  for  work  on 
examining  fluctuations  in  a  waveguide,  as  will  become  apparent 
below.  We  write  the  exact  solution  for  an  unperturbed  case  as 
follows; 


where,  k-  are  the  horizontal  wave  numbers  and  y,-  are  the 
modal  amplitudes  from  which  the  acoustic  pressure  field  is 
obtained.  To  obtain  the  perturbed  wave  functions,  ,  we  seek 
solutions  to  the  equation: 


[fc'(z)-X  Jt/.  =0  , 


which,  for  convenience  we  rewrite  as  follows: 


where 


2=  +  \  .-X.  =  q  +  NX  .  , 


By  integrating  the  overlap  of  the  above  expression  with  'P  j  we 
obtain: 


U  "  u.. 

^  i~  ^u~  .  ‘^ij  where  =  (6) 

<  '  ti 

This  yields  the  Eigenvalue  Correction  Equation 
By  integrating  the  overlap  of  Eq.  5  with  we  obtain; 

j  *  ‘ 

where  jfc  =  1, 2,  3  .  .  .N  and  i  ^  k. 

This  is  an  eigenvalue  equation.  We  can  use  the  expansion 
for  the  eigenvalue  above  to  rewrite  the  equation  as  follows: 

S'  ( -a\  +  H  )  -  I  q  <r..  =  q 

ik\  t  k  tk/  f  1  jk  ij  ik  • 

Where  ik  =  1,  2,  3  .  .  .N  and  i  *  k. 

The  diagonal  terms  are  (a^  ~  +  H  ) 

'  i  k  ik/ 

where  H  q^- AX  .  . 

This  will  prove  useful  later.  The  Hik  terms  contain  the  first  and 
the  higher  order  terms,  some  of  which  are  negligible  in  many 
cases.  The  diagonal  terms  are  almost  always  greater  than  the 
off-diagonal  terms.  Thus  we  may  in  some  cases  use  the  Gauss- 
Seidel  method.  The  first  iteration  via  Gauss-Seidel  leads  to: 


q{z)  =  k^„- k^iz)  and  AX  .  =  X.  -X*  . 

0  III 

For  the  iso- velocity  case  we  can  rewrite  Eq.  1  as; 

‘^"‘*',  2  2  2  0 

- ^  +  a  S'  =0  where  a  =  -  X  .  . 

dz^  '  '  .  0  . 

We  impose  the  following  orthonormaiity  conditions: 

• 

We  assume  closure  so  that  we  can  express  Ui : 
s 

U.=  Ifl.S'. 

‘  y.i 

We  insert  Eq.  4  into  Eq.  3  to  arrive  at: 

Za  ,.  fo*.  .  =  la..  f^(z)  + AX  .V  . 

where i  =  1,2,  .  .  .N. 


“  a^.  -  a^.  +  q..  -  AX  a^.  -a\-  H . 

1  k  ^  kk  I  t  k  li 

N 

where  we  rewrite  ^ ^  ~  ~  ~  • 


This  can  be  a  big  improvement  over  the  ordinary  perturbation 
term  in  that//{jk>  <7iik3nd  thus  we  expect  that  generally  we  will 
have  convergence  if  the  master  matrix  is  diagonally  strong. 
We  see  readily  that  we  can  anive  at  two  improved  perturba¬ 
tion  theories  by  retaining  the  full  form  of  Hik  or  by  excluding 
the  last  term  ( the  higher  order  terms)  of  Hik-  The  advantage 
to  excluding  the  last  term  is  that  we  need  not  employ  a 
sophisticated  perturbation  approach  to  be  discussed  below. 
Further  this  expression  shows  that  the  old  theory  overesti¬ 
mated  the  expansion  coefficient. 

The  new  complete  perturbation  expansion  for  the  eigen¬ 
value  is  now: 

0  N  q  q 

X.  =  X.-q..-  V - 

'  '  “  .  ,.y  a'  -a'.+ W  .  • 
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This  expression  indicates  that  the  first  order  correction  is  still 
the  same  but  that  the  higher  order  corrections  are  over/under 
estimated  in  the  old  theory.  The  preceding  equations  can  be 
reordered  to  form  an  eigenvalue  problem.  The  details  of  the 
solution  of  this  problem  will  be  omitted  here.  Instead,  we 
proceed  directly  to  the  results  in  a  form  most  relevant  for  our 
present  application,  namely  the  expression  for  the  acoustic 
field.  The  acoustic  field  can  be  readily  computed  as  a  sum  of 
appropriate  wavefunctions  obtained  by  the  usual  method 
of  separation  of  variables.  This  leads  to  the  following  expression 
for  the  acoustic  pressure  in  an  acoustic  ductof  water  thickness  h: 


P=pc)  X«/mSina„Z  - 

Im  'Jk  I  r 


ik  I  r  „ 
e  I  .-p,  r 


e~^‘  ^  0<Z<h. 


where  the  aim  and  a„  are  as  defined  previously,  p  is  density, 
(0  is  radian  frequency,  kf  =  is  the  horizontal  wavenum¬ 
ber,  and  P/  is  the  attenuation  coefficient. 


in  the  curve,  appear  to  migrate  vertically  with  time  -  a  motion 
often  suggesting  the  presence  of  internal  waves  (Gregg  [7]). 
Of  particular  note  is  the  fine-structure  evident  in  the  expanded 
portion  of  the  profile  (Fig.  Ic.).  Generally,  the  fine-structure 
in  this  example  was  characterized  by  vertical  dimensions  of 
the  order  of  a  few  centimeters  to  one  or  two  meters.  This  fine- 
structure  tends  to  scatter  high  frequency  sound. 

The  most  salient  characteristics  of  the  sound  speed  profiles 
of  Figure  1  are  the  negative  gradients  of  the  thermoclines.  As 
a  result,  the  sound  fields  will  be  dominated  by  downwardly 
propagating  energy,  resulting  in  substantial  interaction  with 
the  seafloor  and,  most  likely,  channeling  into  a  bottom  sound 
duct.  However,  the  differences  between  the  individual  profiles, 
manifested  by  the  apparent  oscillation  of  the  surface  duct 
thickness,  are  likely  to  lead  to  distinguishable  differences  in 
acoustic  transmission,  even  through  the  general  characteristics 
of  the  propagation  may  be  similar.  It  is  well  known  that  even 
a  relatively  small  change  in  sound  speed  may  lead  to  significant 
alterations  in  underwater  acoustics  propagation,  particularly 
over  long  ranges. 


3.  THE  MEASURED  SOUND  SPEED  PROFILES 


4.  RESULTS  OF  THE  CALCULATIONS 


The  sound  speed  profiles  used  here  were  measured  during 
an  acoustic  propagation  experiment  conducted  in  shallow 
water  (80  m)  in  the  Mediterranean  by  the  S  ACL  ANT  Undersea 
Research  Centre. 

Figure  1  provides  an  example  of  the  variation  of  sound 
speed  with  depth  over  a  peiiod  of  25  hours  (Ali  et  al  [6]). 
Figure  la  plots  in  the  usual  manner  the  sound  speed  profiles 
over  approximately  hour  intervals.  In  Figure  lb  the  same 
profiles  are  displayed  side  by  side  to  clearly  demonstrate  the 
nature  of  the  variation,  while  Figure  1  c  is  shown  to  emphasize 
the  sound  speed  fine-structure.  From  the  surface  down  to  about 
20  m  the  region  appears  to  be  essentially  iso  velocity  with 
sound  speed  of  approximately  1538  m/s.  The  depth  from 
approximately  25  to  35  m  comprises  the  steepest  portion  of  the 
thermocline.  Fluctuations  in  sound  speed  are  quite  evident, 
particularly  at  a  depth  of  about  25  m  within  the  thermocline. 
The  plot  clearly  indicates  an  oscillation  in  the  thickness  of  the 
mixed  layer  (surface  duct).  Further,  “kinks”,  typical  features 


In  order  to  concentrate  on  the  essential  physics  of  the 
propagation  in  the  watercolumn,  we  have  modelled  the  medium 
as  an  80-m  water  column  overlying  an  infinite  half  space 
bottom  with  the  following  geoacoustic  properties: 
compressional  speed  and  attenuation  1650  m/s  and  0.8  dBA, 
respectively,  and  relative  density  1.9.  These  parameters  are 
typical  of  sand-silt  bottoms.  We  have  ignored  shear  effects  in 
the  bottom,  which,  for  the  relatively  high  frequencies  we  are 
considering  here,  play  a  negligible  role  on  the  propagation  in 
the  water  column. 

Calculations  of  the  propagating  acoustic  field — in  particular 
the  mode  shapes  (wave  functions)  and  transmission  loss  as  a 
function  of  receiver  depth  in  the  water  coh'mn  —  were  made 
for  several  different  frequencies  and  source  depths.  To  demon¬ 
strate  the  effect  on  the  acoustic  field  of  differences  in  the  sound 
speed  profiles,  we  could  have  proc'’  Jtied  in  several  ways.  One 
method  is  to  compute  the  appropriate  field  for  each  measured 
profile  and  compare  the  results  for  the  various  profiles.  An 
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alternative  approach  is  to  generate  a  depth-averaged  profile 
from  the  measured  set  and  compare  the  acoustics  of  this  “mean 
profile”  with  that  of  each  of  the  others.  We  chose  this  second 
ai^roach  for  at  least  two  reasons.  The  first  is  mathematical:  a 
key  strength  of  the  mathematical  model  is  the  ease  with  which 
it  handles  a  perturbation  effect.  Therefore,  each  of  the  mea¬ 
sured  profiles  is  considered  a  perturbation  from  the  mean 
profile,  with  readily  computable  coefficients  for  the  normal 
mode  expansion.  Second,  comparison  with  a  mean  proHle  is 
sensible  from  a  physical  point  of  view.  The  mean  profile, 
although  it  does  not  represent  an  actually  measured 
profile,  helps  enhance  our  understanding  of  the  physics  of  the 
propagation.  For,  the  mean  profile  represents  a  “smoothed” 
version  of  the  actual  profiles:  it  eliminates  finestructure  and 
other  “details”  while  retaining  the  essential  large-scale  features 
of  the  propagation  problem.  But,  of  course,  it  is  these  details 
which  distinguish  the  various  profiles;  and  these  effects  are 
more  readily  apparent  if  compared  against  a  mean  profile  than 
against  each  other. 

The  results  shown  here  were  calculated  for  three  profiles: 
the  mean  and  two  perturbed  profiles,  the  latter  chosen  to 
represent  extremes.  The  profiles  are  shown  in  Figure  2; 
Profile  1  is  the  first  profile  displayed  in  Figure  1,  while 
Profile  2  is  profile  28,  almost  at  the  end  of  the  display.  For 
each  of  the  three  profiles  shown,  a  number  of  normal  modes 
was  calculated,  for  a  source  depth  of  70  m  and  a  frequency  of 
3  kHz.  In  Figure  3,  several  of  the  modes  are  compared.  There 
are  several  noteworthy  features.  Independent  of  profile,  the 
lowest  mode  (which  corresponds  to  the  smallest  grazing  angle 
for  the  equivalent  ray)  is  more  concentrated  in  the  region  of 
low  sound  speed  below  the  thermocline,  close  to  the  bottom. 
The  higher  order  modes  (with  larger  grazing  angles)  occupy 
substantially  greater  portions  of  the  water  column .  A  part  from 
the  preceding  behavior,  expected  for  all  the  profiles,  the  modes 
reveal  individual  differences  attributable  to  variations  among 
the  profiles.  It  is  noted  that,  for  the  lower  modes.  Profile  1 
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Figuie  2.  The  dqxh-ivenged  mean  and  two  exlRine  sound  speed  profila. 
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Figure  3.  Selected  modal  functions  at  3  kHz. 
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Figure  4.  Tnuismission  loss  for  a  3-kHz  source  at  a  depth  of 
70  tn  for  (a)  the  mean  profile,  (b)  extreme  profile  I ,  and 
(c)  extreme  profile  2. 


shows  a  greater  deviation  from  the  mean  profile  than  does 
Profile  2.  Though  not  shown  here,  this  behavior  is  consistent 
for  at  least  the  first  6  modes.  For  higher  modes,  the  trend  is  less 
predictable.  Since  the  acoustic  field  is  obtained  in  terms  of  a 
coherent  sum  of  the  modal  functions,  the  behavior  evidenced 
in  Figure  3  suggests  that  Profile  1  will  lead  to  propagation 
characteristics  more  noticeably  different  from  the  mean  pro¬ 
file  than  will  Profile  2.  This  supposition  is  based  on  the  fact 
that  the  lower  order  modes  will  contribute  more  than  the  higher 
modes  to  the  propagation,  particularly  at  greater  ranges,  since 
the  lower  modes  will  be  less  degraded  by  interaction  with  the 
bottom. 

On  the  whole,  the  transmission  loss  contours  shown  in 
Figure  4  support  the  preceding  a  priori  contention.  As  seen, 
channeling  along  the  bottom  is  the  dominant  feature  of  all 
three  figures,  and  the  degree  of  channeling  (i.e.,  the  energy 
levels)  of  the  mean  and  Profile  2  cases  are  very  similar. 
Profile  1  clearly  leads  to  stronger  coupling  with  the  bottom 
channel.  In  retrospect,  one  could  have  expected  this,  based  on 
the  somewhat  larger  negative  gradient  seen  in  Profile  1 
(Figure  2).  An  importantcommon  feature  of  the  two  perturbed 
sound  speed  profiles,  not  shared  with  the  mean  profile,  is  the 
formation  of  a  second  duct  centered  around  45-50  m  depth. 
Clearly,  this  second  channeling  is  attributable  to  the 
quasi-bilinear  feature  evident  in  the  two  profiles  at  approxi¬ 
mately  50  m.  On  the  other  hand,  the  mean  profile  leads  to  a 
fairly  well-defined  surface  duct  centered  around  a  depth  of 
10  m,  a  result  whose  plausibility  follows  from  a  perusal 
of  Figure  2.  A  closer  look  at  the  propagation  behavior  at 
particular  depths  may  be  obtained  by  a  conceptual  horizontal 
slice  of  Figure  4.  The  results  at  two  depths,  45  m  and  75  m,  are 
shown  out  to  a  range  of  30  km  in  Figure  5.  These  figures  show 
even  more  strikingly  the  increasing  divergence  in  transmis¬ 
sion  loss  arising  from  differences  in  the  sound  speed  profiles. 

One  could  continue  in  this  vein,  obtaining  a  good  deal  of 
useful  and  interesting  information  regarding  the  subtleties  of 
sound  propagation  arising  from  variations  in  the  sound  speed. 
But  the  objective  here  was  to  demonstrate  the  utility  of  the 
newly-developed  theoretical  approach  to  handle  problems  of 
this  nature. 


5.  SUMMARY  AND  CONCLUSIONS 

A  newly-developed  high-order  perturbation  method  has 
been  described  and  applied  to  the  problem  of  sound  propagation 
in  a  time- vary  ing  ocean  environment.  Using  a  series  of  measured 
sound  speed  profiles  a  set  of  basis  functions  was  generated  and 
used tocalculate the resultingsound  propagation inashallow- 
water  environment.  The  results  confirm  the  significant  effects 
on  acoustic  transmission  of  seemingly  minor  variations  in 
sound  speed  and,  moreover,  demonstrate  the  efficacy  of  the 
new  perturbation  method  in  handling  such  problems. 
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Figure  S.  Tranimission  loss  at  receiver  depths  of  (a)  7S  m  and  (b)  45  m,  for 
a  3-kHz,  70-m  depth  sound  source. 
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